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MICROWAVE  INTERACTION  WITH  AIR 


I.  INTRODUCTION 

Microwave  breakdown  studies  of  gaseous  elements  have  been  carried 
out  extensively  over  a  wide  range  of  pressures  and  for  several  microwave 
frequencies  using  CW  and  pulsed  radiation  sources1-^  The'  main  emphasis  in 
these  studies  was  on  the  determination  of  the  breakdown  power  threshold  and 
its  dependence  on  the  gas  pressure  and  the  microwave  frequency.  The 

coupling  of  microwave  energy  into  tne  breakdown  plasma  and  neutral  gas  has 
not  been  studied  in  detail.  The  reason  fvr  this  is  that,  until  recently, 
no  high-power  microwave  sources  have  been  available  to  perform  such 
studies.  Most  of  the  early  work  performed  on  breakdown  thresholds  was 

performed  using  high  Q  cavities  to  obtain  the  necessary  electric  field  to 
break  down  the  gas.  Once  breakdown  of  the  gas  occurred,  the  Q  of  the 

cavity  dropped  and  the  interaction  changed.  Using  the  NRL  high-power 

gyrtron  facility,  we  have  been  able  to  eliminate  the  need  for  cavities  and 
have  performed  experiments  using  a  focused  geometry  to  examine  the 
coupling  of  microwave  energy  to  nitrogen  gas  during  breakdown.  We  have 
also  modeled  the  experiments  using  a  1-D  computer  simulation  code^. 
Simulations  were  performed  in  a  spherical  geometry  using  a  self-consistent, 
nitrogen  chemistry,  wave  optics,  microwave  breakdown  simulation  code,  MINI.^ 

The  main  emphasis  of  past  work  was  on  the  ionization  front  created 

during  nitrogen  breakdown  and  its  motion  and  plasma  properties,~a5“observed 

experimentally . -j^This  motion  is  not  a  movement  of  the  plasma,  but  rather  a 

movement  of  where  the  strongest  ionization  of  the  gas  is  taking  place.  We 

believe  it  to  be  due  to  two  effects:  1)  delayed  ionization,  and  2) 

reflection  from  the  plasma  itself.  Since  the  electron  density  is 

tl 

proportional  to  e  °,  where  I0  is  the  local  microwave  intensity, 
breakdown  occurs  faster  in  those  regions  of  high  intensity.  In  a  focused 
system  this  results  in  breakdown  beginning  at  the  focus.  The  breakdown 


then  moves  back  into  the  lower  intensity  region  as  time  progresses.  If 
nothing  else  is  occurring,  one  would  then  ooserve  a  cone  of  plasma  late  in 
time.  However,  the  plasma  is  both  absorptive  and  reflective.  The  reflec¬ 
tion  results  in  standing  nave  patterns.  Thus,  peaks  A/2  apart  are  seen  to 
grow  in  time.  The  first  peak  is  A/4  from  the  plasma  surface.  As  the  dens¬ 
ity  increases  in  these  peaks,  they  become  reflet  ve  and  absorptive, 
cuttina  off  the  microwave  energy  to  peaks  behind  them.  Again,  the  density 
grows  as  eat^°  and,  therefore,  those  peaxs  nearest  the  focus  initially 
grow  fastest.  If  only  one  peak  were  reflecting,  we  would  expect  A/2  struc¬ 
ture.  However,  all  the  peaks  reflect  and,  therefore,  we  exoect  A/4,  the 
distance  from  the  reflecting  surface  to  the  first  electric  field  maximum, 
to  occur  as  well.  Thus,  in  the  experiment  we  would  expect  a  background 
plasma  due  to  delayed  ionization  with  A/2  and  A/4  structures  in  it  due  to 
reflection.-  Reflection  and  absorption  will  also  prevent  microwave  penetra¬ 
tion,  so  we  exoect  to  see  the  plasma  disappear  behind  the  front  of  the 
ionization.  This  results  in  the  appearance  of  motion.  This  motion  may 
then  occur  as  A/2  or  A/4  jumps,  since  once  tnese  structures  develop,  the 
intensity  at  these  points  is  larger  due  to  reflection  than  the  background 
intensity. 


This  ionization  front  motion  strongly  affects  the  coupling  of 
microwave  energy  to  the  gas;  it  determines  where  ti*e  energy  will  be  ab¬ 
sorbed,  since  its  nigh  density  does  not  allow  microwave  penetration.  Also, 
since  it  is  moving,  it  spreads  the  heating  out  in  the  gas.  Motion  towards 
the  microwave  source  has  been  observed  by  Salisoury  and  Fiint4,  and  by 
Beust  and  Ford5.  The  velocities  were  much  slower,  however,,  V  *  3  x  103 
cm/s-  for  Salisbury  and  Flint's  work,'  and  V  «  600  cm/s  for,  Beust  and  Ford's 
work.  Raizer6  advanced  a  theoretical  treatment  for  the  microwave  discharge 
propagation  in  high  pressure  air  (=1Atmj  to  explain  Beust  and  Ford's  re¬ 
sults  under  CW  conditions.  Other  work  by  Scharfman,  et  al7,  observed  that 
for  low  pressures  the  breakdown  plasma  became  opaque  to-  microwaves,  thus 
preventing  penetration.  In  the  previous  study,,  we  observed  the  formation 
of  the  front,  its  motion,  the  prevention  of  microwave  penetration,  and 
measured  the  plasma  density  and  electron  temperature  of  the'  front.  Good 
agreement  was  found  between  the  experiment  and  the  simulation  model.  This 
work  is  briefly  reviewed  in  Section  II. 
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II. 


EXPERIMENT 


A.  Purpose  of  the  Experiment 

MRC  has  participated  in  establishing  an  experimental  air  breakdown 
facility  at  NRL.  In  the  past,  a  35  GHz,  150  kV  gyrotron  was  used.  A  200 
pps  repetition  rate  and  2-  ys  pulse  length  were  possible.  Microwave  pulse 
shape  and  power  were  very  stable.  Breakdown  effects  also  seemed  to  be 
stable,  thus  allowing  data  to  be  taken  over  many  pulses.  The  microwave 
power  available,  however ,. limited  breakdown  to  a  pressure  of  approximately 
100  Torr.  A  newly  developed  95  GHz,  1  MW,  13  ys  pulse  length  gyrotron  was 
to  be  made  available  ? .  a  microwave  source  for  this  facility  and  would 
allow  breakdown  of  air  at  atmospheric  pressure.  The  air  breakdown  program 
was  to  be  completed  by  performing  experiments  at  full  atmospheric  pres¬ 
sure.  This  would  allow  one  to  determine  if  any  physical  processes  become 
important  due  to  increased  collision  frequency.  Further,  the  atmospheric 
pressure  region  is  the  one  for  which  microwave  weapons  have  been  con¬ 
sidered. 


The  experimental  plan  called  for  experiments  to  be  performed  to 
determine  the  microwave  coupling  to  air  at  one  atmosphere.  In  particular, 
emphasis  was  to  be  placed  on  determining  if  the  ionization  front  behavior 
continues  and,  if  so,  how  it  varied  with  pressure.  Further,  a  determina¬ 
tion  to  see  if  the  ionization  front  motion  could  be  prevented  was  to  be 
made.  In  order  to  observe  the  ionization  fronts,  high  speed  photography 
was  planned.  The  plasma  density  and  temperature  were  to  be  monitored. 
This  was  to  be  done  using  existing  diagnostics  developed  by  MRC  for  the  NRL 
air  breakdown  facility  and  modified  to  improve  their  performance. 

Since  the  planned  research  was  an  extension  of  previous  work,  the 
previous  experiments  and  results  are  reviewed  below. 

B.  Previous  Experiments 

The  microwave  energy  coupling  experiments  were  performed  in  dry 
nitrogen  using '  one  of  the  HRl  35  GHz,  high  power  gyrotons8.  This  gyrotron 


is  capable  of  150  kW,  lysec  pulses  at  100  pps.  The  investigation,  however, 
was  performed  with  a  microwave  power  =  112  kW  using  S  polarization.  The 
nitrogen  pressure  was  variable.  The  data  presented  is  at  a  pressure  o^  25 
Torr  which  corresponds  to  v/w  -  0.65,  where  v  is  the  collision  frequency 
for  electron  momentum  transfer,  and  u  is  the  natural  frequency  of  the 
microwave  radiation.  The  microwave  radiation  is  introduced  into  a  test 
chamber  using  a  conical  horn  with  a  focusing  dielectric  lens  in  front  (see 
Figure  1).  The  lens  has  a  focal  length  of  11.2  cm  and  a  diameter  of 
7.62  cm.  The  focal  spot  is  elliptical  and  has  an  approximate  Gaussian  dis¬ 
tribution  in  intensity,  with  half  widths  of  1.5  cm  and  1.1  cm.  The  3dB. 
focal  spot  area  is  approx:  nately  1.7  cm2.  For  112  kW,  this  implies  33 

o 

kW/cm  average  power  inside  the  3dB  spot.  Breakdown,  without  a  reflecting 
surface,  was  observed  to  cease  at  p  =  75  Torr  fee  P  =  112  kW.  Theoretic¬ 
ally,  we  calculate  an  average  power  of  27  kW/cm2  is  required  for  breakdown 
at  75  Torr  in  nitrogen.  This  is  in  reasonable  agreement  with  an  estimated 
intensity  of  33  kW/cm*.  We  have,  however,  ignored  the  gaussian  spread  and 
assumed  an  even  distribution. 

The  breakdown  experiments  were  performed  with  and  without  a  metal 
surface  at  the  focal  point.  The  planar  met?.;  surface  couid  be  oriented 
perpendicular,  90°,  or  with  an  angle  of  45*  to  the  incident  radiation.  The 
metal  surface  serves  to  set  uy  a  standing  wave  pattern  in  front  of  the 
surface,  resulting  in  4  I0  at  the  peaks,  where  I0  is  the  intensity  of 
the  incident  radiation.  Also,  a  thin  dielectric  could  be  placed  off  the 
metallic  surface.  This  was  used  in  determining  the  cause  of  the  plasma 
motion  observed,  and  will  be  discussed  later. 

Observation  of  the  plasma  was  done  using  a  TRW  framing'/ streaking 
camera.  A  20  ns  exposure  time  was  used  in  the,  framing  mode,  aind  a  2  ys 
exposure  time  in  the  streaking  mode.  Due  to  the  amount  of  light  present, 
multiple  exposures  were  required  to  obtain  easily  visible  photographs. 
However,  some  single  shot  experiments  were  performed  and  showed  that  the 
results  were  reproducible  shot-to-shot.  The  electron  temperature  and 
density  were  measured  by  using  a  nitrogen-helium  gas  mixture.  In  a  pure 
nitrogen  mixture,  only  second  positive  nitrogen  bands  were  observable.  By 
introducing  (ns^ium,  the  5876A  line  could  be  used  with  the  nitrogen  337 1 A 
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Figure  1 .  Schematic  of  nitrogen  breakdown  experiment 


band  to  make  the  required  density  and  temperature  measurements.  Ar,  Oriel 
model  7240  monochrometer *  f/3.7,  with  a  10A  bandwidth  was  used.  The  mono¬ 
chrometer  output  was  measured  using  a  102°.  photomultiplier  with  a  base  mod¬ 
ified  for  nano-second  response  time.  Absolute  intensity  calibrations  were 
performed  using  a  calibrated  deuterium  arc  lamp  source9  for  3371 A  and  a 
calibrated  strip  tungsten  lamp  source10  for  5876 A. 


Our  experiments  were  performed  at  25  Torr  where  we  had  approxi¬ 
mately  nine  times  the  intensity  required  to  brea.K  down  the  nitrogen  at  the 
focal  spot.  With  the  presence  of  a  metal  surface,  a  factor  of  four  reduc¬ 
tion  in  the  power  required  for  breakdown  occurs  due  to  the  standing  wave 
pattern  which  is  set  up.  Due  to  the  focused  nature  of  the  experiment,  the 
standing  wave  pattern's  intensity  drops  as  the  distance  from  the  surface  is 
increased.  For  the  45°  target,  in  fact,  there  is  only  a  small  region  of 
standing  waves  perpendicular  to  the  surface.  It  ^s  only  for  the  90°  sur¬ 
face,  or  the  case  of  a  reflective  plasma  perpendicualr  to  the  microwav*- 
beam,  that  a  standing  wave  pattern  back  the  length  of  the  beam  line 
exists.  A  delayed  breakdown  from  the  surface  back  towards  the  lens,  is  then 
expected,  since  the  growth  of  the  plasma  is  proportional  to  e^0.  We 


believe  that  fluctuations  in  the  reflective  plasma  surface  and  UV  ioniza¬ 
tion  (if  any)  will  cause  a  smearing  out.  of  the  standing  wave  pattern. 


1 .  Time  Resolved  Photographs  of  Breakdown 

The  observed  breakdown  patterns  fbr  the  three  cases  are  shown  in 
Figures  2  througn  4.  In  rigure  2,  the  breakdown  without  a  surface  is 
shown.  A  conical  plasma  is  formed  with  its  tail  nearest  the  focal  point 


and'  head  towards  the  lens..  This  indicates  that  a  large  fraction  of  the 
.microwave  energy  is  concentrated  on  axis.  For  a  diffraction  limited  spot 
size,  we  calculate  the  depth  of  field  (where  I  *  1 0/e)  for  our  case  to  be 
3.5  cm,  which  is  approximately  the  length  of  the  tail  in  the  early  break¬ 
down  region  observed  in  Figure  2.  The  total  length  of  the  breakdown  region 
is  5  cm.  We  calculate  a  length  of  3.7  um  based  oft  incident  intensity 
decreasing  to  the  nitrogen  breakdown  threshold  as  we  move  away  from  focus. 
This  suggests  reflection  off  the  plasma  must  be  occurring.  We  caicualte 

that  *  19%  reflection  would  be  required;  the  simulation  suggests  R'  >_  29%. 

\  .  •  '  '  * 

6  .  .  ' 


Figure  2.  Framing  photographs  of  breakdown  at  25  Torr  N2 
for  no  surface*  presen c,  I0  x  112  kW  or  33  k\V/cm2  at  focus. 

The  microwaves  are  incident  from  the  left.  Frames  are  100  ns 
apart  and  have  a  20  ns  exposure  time.  Time  shown  on  right 
refers'  to  the  delay  after  the  first  observable  breakdown  and 
not  the  microwave  pulse. 
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Figu-e  3a.  Framing  photographs  at  25  Torr  N2  for  90* 
metal  surface  present  at  focus,  I0  =  112  kW.  Times  on 
right  refer  to  the  del'  after  first  oDservabl^  breakdown. 


8 


r  Igurc  3b.  Streaking  photograph.  at  ?5  Torr  ll2  for  90°  metal  surface 
present  at  focus,  I0  =  11?  kW.  Streaking  time  was  2  us.  The  velocity  of 
motion  for  the  ionization  front  Is  measured  to  be  4.6  x  1 0 6  cm/s  from  the 
photograph.  - 


% 


Figure  4.  Framing  photographs  at  25  Torr  M2-  for  45* 
metal  surface,  I0  =  112  k*.  Times  shown  ort  Fight 
refer  to  the  delay  after  first  observable  breakdown. 
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The  breakdown  begins  at  tr>c  focal  spot,  but  since  the  microwave  power  is 
above  the  breakdown  threshold,  a  secondary  breakdown  region  ahead  of  the 
focus  forms  in  a  time  *  2.4  x  10*  ;s  after  the  initial  breakdown.  As  the 
electron  density  in  this  secondary  breakdown  region  increases  to  Nc  =  (u* 
*  v  )me/(4i>e  ),  the  tail  (focus)  reqion  disappears  due  to  the  microwaves 
being  absorbed  and  reflected  by  the  plasma  at  the  head.  Cnee  the  microwave 
energy  is  stopped  from  reaching  the  plasma  in  the  tail,  the  plasma  temper¬ 
ature  there  drops  rapidly  (t  «  10  ns),  and  the  emitted  light  ceases.  The 
plasma  itself  also  begins  to  decay,  but  on  a  longer  time  scale.  Late  in 
time,  we  see  that  a  single  absorbing  layer  about  one  wavelength  tnick  is 
formed.  ' 


Figure  3a  shows  the  framing  camera  data  for  the  perpendicular 
surface  experiment,  and  Figure  3b  shows  a  streak  photograph  for  the  same 
case.  We  expect  a  standing  wave  to  be  set  up  off  the  surface  with  the 
first  maximum  X/4  from  the  surface,  and  with  the  rest  of  the  maxima  spaced 
X 12  apart.  This  can  be  seen  in  Fiqure  3a.  Also  observed  is  a  "blooming" 
of  the  visible  plasma  near  the  tarqet.  We  believe  this  is  an  effect  due  to 
the  focusing.  The  microwave  beam  becomes  more  intense  near  focus.  There¬ 
fore,  the  diameter  of  the  beer,  capable  of  breaking  down  toe  plasma  in¬ 
creases.  .  Similar  to  the  no-target  case,  an  absorbing  front  travels  away 
from  tne  tarqet  towarus  the  microwave  source.  From  the  displacement  in  the 
streak  photograph,  we  calculate  that  the  plasma  moves  away  from  the  surface 
with  a  velocir>  »  4.6  x  10fe  cm/s.  It  continues  to  move  back  along  the 
microwave  beam  until  the  beam  energy  decreases  sufficiently  (due  to|de- 
focusinq)  so  that  the  beatinq  with  any  reflected  microwaves  is  below  the 
breakdown  threshold.  This  "stagnation  point"  is  also  observed , for  the  45* 
surface.  Similar  behavior  is  also  observed  in  the  computer  simulation. 

The  ^5*  data  is  snown  in  Figure  4.  It  is  similar  to  the  90*  data, 
except  fringes  occur  X/(?  cos  6)  apart,  b  =  i5*,  with  the  first  fringe  at 
X / ( 4  cos  -)'  from  trie  surface.  Again,' an  absorbing  front  moving  alongj  the 

beam  axis  is  observed.  This  occurs  as  a  result  of  enerc ■  reflected  by|  the 

plasma  along  tr>e  beam.  Since  we  are  far  above  threshold,'  e  belieye  a 

planar  plasma,  perpend icular  to  the  beam,  qrows  near  the  target  anrjl  i: 
initiated  by  the  plasma  growth  at  the  standing  wave  maxima.  Its  density 


yrows  (Ng- ■— w— Nc) ,  and  it  becomes  reflective,  setting  up  standing  waves. 
Plasma  density  builcs  at  these  maxima,  cutting  off  microwave  penetration, 
and  the  process  repeats.  Visually,  this  appears  as  plasma  motion.  be 
would  expect  the  motion  to  jump  in  either  A/2  or  A/4  steps.  Some  of  the 
45°  data  shows  A/2  bars  behind  the  main  moving  plasma.  In  general,  how¬ 
ever,  this  is  not  observed.  The  A/2  spacing  is  at  about  the  limit  of  reso¬ 
lution.  Further,  we  believe  that  fluctuations  in  the  reflective  plasma 
surface  cause  a  smearing  out  of  the  standing  wave  pattern  and  results  in 
the  more  solid-like  appearance  we  observe  for  the  moving  plasma. 

One  other  possible  explanation  for  the'  front  motion  is  that  U 
emission  causes  ionization  to  occur  in  front  of  the  plasma.  This  new 
plasma  absorbs  microwaves,  heats  up,  and  then  the  process  repeats.  This 
would  result  in  a  uniform  smear  with  no  A/2  structure.  Our  data  is  some¬ 
what  suggestive  of  this.  To  check  this  idea,  we  placed  a  thin  dielectric  a 
few  wavelengths  of'  the  metallic  surface  to  stop  UV  penetration,  but  allow 
microwave  transmission.  If  UV  ionization  were  the  method  of  propagation, 
we  expected  the  plasma  to  stop  moving  when  it  hit  the  dielectric.  The 
plasma  moved  through  the  dielectric  barrier.  be,  therefore,  believe  the 
motion  to  be  due  to  reflection  from  the  plasma  itself,  resulting  in  ioniza¬ 
tion  fronts.  However,  ionization  due  to  UV  may  be  responsible  for  some  of 
the  smearing  of  the  A/2  and  A/4  structure.  The  UV  ionization  mean  free 
path  is  -  0.05  cm  «  A/2  for  our  conditions. 

For  the  45*  target  case,  t!it  breakdown  was  examined  as  the  power 
was  decreased.  As  the  power  was  decreased,,  the  ionization  front  was  ob¬ 
served  to  staqnate  closer  and  closer  to,  the  target,  as  expected.  When  the 
incoming  power  was  20  kW,  •  5.9  kW/cm  ,  we  observed  only  two  fringes  and  no 
apparent  ionization  front.  At  30  kW,  the  front  reappeared.  Incident  radi¬ 
ation  on  a  metal  surface  for  breakdown  to  occur  at  25  Torr  is  -  0.75 
, kW/cm2.  However,  the  pressure  was  only  twice  the  free  nitrogen  breakdown 
threshold  of  3  kW/cm^,  arid  this  will  be  tne  dominant  threshold  for  ioniza¬ 
tion  front  travel  due  to  delayed  ionization. 


2 .  Spectroscopic  Measurement  of  Density  and  Temperature 


"n  order  to  measure  the  temperature  and  density  of  the  electrons, 
the  absolute  intensity  of  the  nitrogen  second  positive  band,  3371  A,  was 
compared  to  the  absolute  intensity  of  the  neutral  helium  line  at  5S76A.  A 
mixture  of  the  helium  and  nitrogen  was  used  to  perform  this  experiment. 
The  3371A  fcano  corresponds  to  the  (0,01  transition  in  the  second  positive 
band  system  of  N,.  ^e  obtained  the  electron  temperature  from  the  ratio  of 
the  two  lines1 1  and  the  electron  density  from  the  absolute  intensity 
measurement  and  the  deduced  electron  temperature. 

For  an  optically  thin  line,  and  assuming  a  steady  state  coronal 
model,  the  band  intensity  at  3371 A  can'  be  expressed  as: 

N  N-  X  A  E  L 

..  e  2  oc  o,o  o,o  ,  Ml 

Uui  - - - - * - 4  '-q  ~7T- 

o  2 

where  Nt  is  the  density  of  the  nitrogen  molecules  in  the  ground  state, 
\oc.  is  the  emission  excitation  rate  coefficient,  ,A00  is  the  transition 
rate  for  the  band  at  3371 A  and  A0  is  the  total  transition  rate  for  the 
excited  state,  E0>0  is  the  pnoton  energy,  L  is  the  length  of  the  emitting 
plasma,  and  q  is  the  quenching  rate  coefficient  of  the  excited  state  by  the 
neutral  species.  A  similar  equation  can  be  written  for  f.a  intensity  of 
the  helium  line.  The  excitation  rate  coefficients  for  3371A  and  5S76A  were 
obtained  using  the  measured  cross  sections12*1^  averaged  over  a 
“Maxwellian  electron  velocity  distribution.  The  transition  rates  for  3371« 
and  for  5S7oA  are  from  Refs.  (14)  and  (15),  respectively..  The  quenching 

rate  coefficient*6  for  3371  A  by  N2  is  1.15  x  10-  11  cm5/sec. 

The  measured  absolute  intensities  of  3371A  and  5576A  in  a  mixture 
of  15  Torr  M2  and  19  Torr  He  were  25.36  %/cm2sr  and  5.6  x  10-2  W/cm2sr. 

Using  the  above  analysis,  we  obtain  Te  =  3.5  eV  and  Me  =  4.4  x  10 

cm”^  for  the  average  density  over  a  2  cm  chord.  If  quenching  of  the  helium 
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line  by  nitrogen  is  included,  the  electron  temperature  changes  slightly, 
Te  =  4.1  eV,  with  a  corresponding  slight  decrease  in  the  calculated  elec¬ 
tron  density.  This  electron  density  is  in  reasonable  agreement  with 
Langmuir  probe  data,  Ne  "  2  x  10i3  cm-3,  and  simulation  predictions,  Ne 
=  3.6  x  10i3  cm-3,  since  these  last  two  are  peak  local  densities.  It  is 
reasonable  to  expect  that  the  peak  density  is  at  least  four  times  larger 
than  the  average  density,  since  the  density  is  proportional  to  e  ot, 
and  the  microwave  intensity  is  Gaussian  across  the  chord.  (The  two  centi¬ 
meter  chord  is  slightly  larger  than  the  FWHM  of  the  intensity  Gaussian.)  A 
comparison  of  the  absolute  intensity  of  the  337! A  nitrogen  band  for . pure 
fJ2 ,  as  observed  experimentally  (Iexp  =  1*12  x  1020  eV/ (cm2-s-sr ) ) ,  agreed 
reasonably  well  with  the  simulation  predictions  for  similar  conditions' 

o  n  o 

( I ,im  =  3.6  x  10  eV/(s-cm  -sr)).  As  expected,  the  simulation  also  pre¬ 
dicted  a  lower  electron  temperature,  Te  »  3.25  eV,  for  the  pure  N2  case. 

In  conclusion,  we  have  observed  that  for  I  >  lg,  ionization 
fronts  form  and  move  towards  the  source  at  a  velocity  -  4.6  x  196  cm/s. 
These  ionization  fronts  result  from  both  delayed  ionization  and  reflection 
of  microwave  energy  from  the  front.  The  front  density  can  become  quite 
large  and  is  seen  to  cut  off  microwave  penetration  in  about  one  wavelength 
in  late  time.  The  average  chord  piasma  density  was  spectroscopically  mea¬ 
sured  to  be  Ue  =  4.4  X‘1012  cm-3,  and  the  electron  temperature  to  be  Te 
=  3.5  eV  for  a  nitrogen-helium  mixture.  Langmuir  probe  measurements  of  the 
density  gave  reasonable  agreement  with  peak  densities  for  pure  nitrogen  of 
Ne  ■  2  x  10 13  cm-3.  The  cut-off  of  microwave  penetration  observed  photo¬ 
graphically  suqgests  an  electron  density  Ne  =  Nc  =  2.2  x  10l3  cm"3.  Vie 
expect  the  electron  temperature  in  the  gas  mixture  to  be  higher  since  N2 
vibrational  and  low-lying  electronic  states  act  to  decrease  Te,  and  mono- 
acomic  He  has  no  such  states.  Reasonable  agreement,  was  seen  „ith  the 
simulation  for  the  velocity  of  the  ionization  front,  Ne,  Te,  and  I337i, 
The  simulation  predicts  small  reflection,  R  <  29%,  and  low  gas  temperature, 
Tg  <  0.027  eV. 
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C.  Present  Experiment 


Several  modifications  to  the  previous  experimental  set-up  were 
planned  to  improve  results.  In  particular,  the  photomultiplier  tube  used 
in  the  previous  study  was  limited  in  its  response  to  the  5576A  He  line.  A 
new  tube  was  proposed  to  NRL  for  use  to  enhance  the  diagnostic  sensitiv¬ 
ity.  The  same  streaking ,  and  framing  set-up  was  planned  for  use  and  was 
recalibrated. 

Unfortunately,  a  larqe  difficulty  arose  after  the„  plasma  diagnos¬ 
tics  were  configured.  The  IGW  gyrotron  tube  was  not  available  for  use  due 
to  dif f iculties  in  getting  the  tube  operational,  and  WC  began  to  assist  in 
attempting  to  get  the  tube  operational  by  participating  in  the  gun  design. 
These  efforts  are  described  in  the  next  section.  .To  date,  the  IGW  gyrotron 
has  not  been  functional.  All  efforts  to  get  the  device  functional  at  IGW 
have  failed.  , 


III.  SIMULATION  OF  THE  GUN 

A.  Introduction 

This  section  describes  the  results  of  computer  calculations  of 
electron  trajectories  performed  in  support  of  the  gyrotron  design  effort  at 
URL.  The  goal  of  these  calculations  was  to  determine  acceptable  combina¬ 
tions  of  magnetic  field  profile,  intermediate  anode  voltage,  and  anode 
voltage  to  place  the  electrons  at  the  proper  guiding  center  radius  and  with 
the  proper  transverse/axial  velocity  ratio  (alpha)  to  maximize  the  gyrotron 
Interaction  with  the  TE^  waveguide  mode.  The  code  used  in  these  calcula¬ 
tions  is  SCRIBE* 7 ,  a  version  of  the  SLAC  electron  optics  code  originally 
writteh  by  W  Herrmansfeldt.  The  modifications  in  this  version  are  due  to 
R..H.  Dackson  of  MRC.  , 


B.  Simulation  Parameters 

The  simulations  were  performed  within  the  following  parameter 

context: 


maximum  anode  voltage: 

70  kV 

maximum  interelectrode 

voltage: 

35kV 

beam  current: 

10A 

desired  alpha: 

1.5 

drift  tube  radius: 

0 .48cm 

minimum  magnetic  field 

at  cathode: 

2.3kG 

maximun  coil  current: 

50A 

The  major  parameter  to  be  determined  is  the  average  electron 

guiding  center  radius.  Since  the  gyrotron  is  designed  to  interact  with  the 

TEj3  waveguide  modes,  the  electrons  must  be  placed  at  the  proper  radius  for 

maximum  interaction  with  this  mode.  It  has  been  shown18  that  the  inter- 

action  strength  is  proportional  to  the  quantity  (3n_-|  (x)  +  3n+->  (x) ) , 

where  3m  is  the  Bessel  function  of  the  first  kind  of  order  m  and 

x  =  k^r/rf);  k-|  3  =  third  root  of  the  derivative  of  the  first-order 

2  2 

Bessel  function,  and  rO  is  the  drift  tube  radius.  A  graph  of  3q+32  ^or 
3<r<.5cm  is  shown  in  Figure  5.  This  graph  shows  .that  the  maximum  inter¬ 
action  strength  between  the  TE^  mode  and  a  hollow  electron  beam  is  at  a 
radius  of  approximately  0.19cm.  This  is  the  desired  average  guiding  center 
radius.  Adiabatic  considerations  imply  that  a  magnetic  field  compression 
ratio  of  approximately  12:1  (final=33..6kG:cathode=2.8kG)  ,  will  place  the 
electrons  near  the  desired' radius. 


c. 


Simulation  Input  Data 


Boundary  data  for  SCRIBE  input  were  generated  from  a  digitization 
of  the  electrode  configuration  as  specified  in  the  electron  gun  blueprint 
supplied  to  MRC  by  NRL.  For  the  purposes  of  this  simulation,  the  straight 
section  of  the  anode  was  extended  to  permit  trajectory  calculations  in  the 
uniform  magnetic  field  region.  No  other  boundary  or  electrode  shape 
changes  were  made.  Please  refer  to  the  typical  trajectory  plots  shown 
below  for  details  of  the  simulation  boundary. 

The  mesh  sizes  selected  for  this  simulation  were  1/2  and 
1/4  mm/cell:  Computation  speed  was,  of  course,  faster  with  the  coarser 
mesh,  with  results  very  ciose  to  those  obtained  with  the  fine  mesh.  The 
majority  of  the  graphical  results  presented  here  were  obtained  using  tne 
1/4mm  mesh.  A  listing  of  the  boundary  input  data  for  the  1/4mm  mesh  in  the 
standard  SLAC  code  format  is  given  in  Table  1. 


Magnetic  field  data  for  this  simulation  were  generated  using  a 
separate  magnetic  field  code  (EFFI1*3)  that  directly  calculates  the  vector 
potential  at  any  point  in  space  due  to  an  assembly  of  finite  dimension 
coils.  This  approach  is  superior  to  the  conventional  method  of  first  cal¬ 


culating  the  on-axis  fie,ld  due  to  an  assembly  of  coils,  and  then  using  an 
expansion  to  calculate  tie  field  off-axis.  The  manufacturers'  coil  dimen¬ 
sion  and  winding  data  were  used  in  these  calculations.  A  diagram  of  the 
individual  coil  position:;  in  trie  same  coordinate  system  used  in  the  simula¬ 
tion  is  shown  in  Figure  >,  with  the  trim  coil  identified.  The  approximate 
location  of  the  cathode  emitting  surface  is  also  indicated.  .  A  graph  of  a 
typical  calculated  axial  magnetic  field  oro^ile  on-axis  for  the'  main  coil 
assembly  is  shown  in  Figure  7,  in  this  case  driven  at  45. 41  A.  Also  shown 
is  the  field  profile  (multiplied  by  -10)  of  the  trim  coil  driven  at  24.22A. 
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0,  53,  *0,  0.2972336,  P.ooooooo  2,  33. 1  "3,,  (,.7022519 ,  2. 

«—53— 9.1— D-J  20  9X38— X.  6811 98 1 ?— 37j_l  8jt,_!L.  o**in7  j,_2^ 

9,  52,  92,  C. 9036378,  2.0000000  2,  33,1*5,  0.1»5*627,  0. 

o,  52,  93,  0.766*3*1,  2. 0000*00  2,  32,1*6,  0.927*1*5,  ?. 

. 9-— S2— J 9— X59_0Jl3iX_X nnonnon  2.  3?.l»7.  (i.66987ni ,  2. 
«,  52,  95,  0.01323*5,  2.00oo"00  2,  32,1*8,  0.9115257,  2. 

9,  52,  06,  0.2560388,  2.0000000  2,  32,189,  0.1533775,  0. 

9  , — 52-— OX, — 0-US36028-— D-3373813 2«— 0,1*0.  q.8»s2352.  2, 


u  , 

51,  98, 

0. **2*950, 

?.00.00A|>0 

2, 

0.63708*0 

a  f 

51,  99, 

o'.  7060932, 

p.oaooooo 

2. 

3t,t9?r 

o.a?b»3«%, 

_5X_SA— X52A2935,. 

-2,0040006 _ 

2,. 

3l,-<>3.  (,.35*7532 

SI,  51, 

0.3S?aa?5, 

?. 0000090 

2, 

31,190, 

ft.290K0*b, 

51,  5?, 

0.1 ^5*477 , 

0.*93«*«S 

31.1*5, 

G«i2>257*0  ( 

<U- 

-  SfU— S3*. 

-0. 99AAaA<».*. 

?.OO#990fl_ 

2,. 

a. 

5*,  Sa, 

?. Oconto* 

2. 

31,1*7, 

ft .  o  a  ft  a  0  ©  m  , 

«* 

so,  r;f 

0.*a5?a»a, 

?.9009©f»9 

?. 

31, i»*. 

*3i  •?, 

4*. 

0000090 — _ 

2,. 

30,1 8*. 

a , 

SO,  *7, 

0.?^1 *S?7, 

>.009090* 

y. 

JO, .268, 

0.8*21 **9, 

S«, 

ft. 1 1 «#5>9, 

0.«*0  96?77 

2, 

3»,2»1 , 

0.*1 C05*«( 

*  . 

.  ©o. 

ft  ft  ft  GAf-ft 

2. 

30,292, 

,o.T«5»5*0( 

0000000 

nooo  no  a 
nooo«00 

o  o  o  0  ft  ft  ft 
Tpsn*'* 
OOOOOOC 
n  .1  o  o  ft  n  n 
nnooftoo 
OOOOftOO 
o ft on ft ftft 
ftOftflftftft 

oooooor 

fiBflOflftft 

baanono 

6*aoft0ft 

Shannon 

SAtlftftftft 

OC'OOOOG 
0000000 
nnnonno 
0000000  _ 
S  200000 
52nooon 
5200000  . 
S?ft0ft0.1 
52«0ft00 
55253b0_ 
7aS9Sft« 
1313171 
7  5575?#*_ 
73*37?* 
oonoono 
Ol»00-'00_ 
OOftOftOft 
0ooo?o5 
AOftOftOi?— 
ftOftftftOO 
07?«nnA 
ooooooo_ 
ftftftftftftO 
00*0000 

*u6«**hi _ 

OftOOftftO 

OOOOftOO 

00000ftft_ 

720367* 

000000* 

oonoooo 

OOOOftOO 

59fl1b?o 

OOOOCOO._ 

ftOOOft'nft 

O0OOO00 

0000000 _ 

0000000 

noooo*0 

OOOOftOO _ 

ftftftftftftft 

>^90030 
noo^on _ 

OOftOfUtO 

Ofk*AO0* 

OOOOOOO 


Tabie  1,  Scribe  Boundary  Input  Data  For  Gyrotron  Gun 
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2. 

3n,?n3, 

0.6778031 , 

?.  OftftO^OC 

2,  ?*,2**,  0.3*ft«ft*ft,  2.(l(M'0.7'in 

2. 

22,320, 

0.0«893?3, 

2. ooooono 

2, 

30 , 904 , 

ft. 60O7QA5, 

?,Aonoftoo 

p,  2#3,?^7,  0.3‘>03?Ol,  P.0000OAO 

2, 

22,330, 

0.030*437 , 

P.«53oopo 

-2a— 30. -gflSu.- C-.-5(l-l?tl09._2.  0000000 — 

2,  2b,2*ft.  L'^  2522335,  2.(10000(10  _ 

-2.. 

-21.311, 

__0.  9*27-571, 

-2*0000000 

2, 

30 , P06 , 

ft.«73*273. 

P.noonftoo 

?,  2ft, ?*»,  0.1ft414*ft,  ?  f ******* 

?; 

?1  ,^3P, 

n.  *9a^7f.5. 

2  #  oooonon 

?, 

30, POT, 

ft.  4055307, 

P.ftftTOOOft 

2,  2* ,  27  ft ,  (1.  ,*0**3,  P.miftonon 

2, 

21,333, 

0.6P6S?39r 

C . 00  TftOftO 

-i.-JS'USlLP.,. 

_ft— 33744*1_,_ 

_2*008060C— 

2,  2*,  2 71.  .0. (*4707ftfa,  0.7*45*4* 

2. 

21.334. 

—0.3  58  497-2.— 2.  ft  COOOoO 

2, 

30,?oo, 

ft  .  ?ftft3ft  1  *5, 

?,  aoooooo 

?,  ?5 , ?7  2 ,  O.970ftp5l,  p, nooooon 

?, 

21.335, 

0.69nai ftfe. 

2. OOCO^OP 

*' . 

50,21ft. 

0.3012709. 

?.  0060000 

2,  25,?73,  0.41174*5,  2.0000000 

?, 

?•  1^36, 

O.*223?40, 

2#onooopo 

-2— 

—3  0.  21 1_^ 

ft , 1331*8?.  2-Q0Q0Q0ft 

2.-25,274.-0.8437  1  1ft.  2.0000000 

2.  _21.-337. 

-0.5542355. 

2, ftOOOOOO 

2, 

30,21?, 

0.065101b, 

0.45*14*2 

2,  c. 7756292,  2.0000O0O 

?. 

?! 

o.«*61 0*6, 

P.onnmnn 

?. 

2  ,213, 

0.9470190, 

2.0000000 

2,  25,276,  ft.70753».o,  ?.O0ftftOftO 

2, 

2l  ,774, 

ft.ftl 806P2, 

2  .  o  ft  o  o  o  o  o 

-2-.- 

-25*^10*. 

A2B  ft  2h  5  . -JL.  0  0  0  0  D  0  0 

_  2*  -25,27  7,_i^b3«',5«rD,  2*0000000 

2. 

21*  Jan* 

-  Ci.3ft.ft°75fa. 

.2.0000000 

?. 

?® ,  ?  1  5 , 

0. 860830R, 

2.0000000 

?,  ?5 , 27  * ,  0 . 57 l 3  - 3« ,  2,0000000 

P  # 

1  , 

0.2*1 **®U, 

2,OOOftOOO 

2. 

>9, PJb, 

(1.7027532, 

2,0000000 

?,  ?^#?70,  n,503'*65,  ?.ftftftftftftO 

?, 

21 , 3a?  t 

ft . PKHOoon  , 

p#  ononooo 

-2.- 

-2.9*  24  7.  ft»Jl?4bbbb*-2.  OOnOOOO  — 

2  25.  28  0. .  11.4  351  Or*,  2.0000000 

?. 

21.  343, 

.  L*268D0P(j-« 

2.0000000 

2. 

?4,?1», 

n.bftftSftftft, 

2. nooooop 

2,  ?S,?H1,  ll.36>i03o#  2,0000000 

?, 

, «i r , 

fl.?6  9nftft(i. 

p.ftftoftooo 

?. 

?P,21Q, 

0.5*80933, 

2.0000000 

2,  ?5,2H2,  ft.?9'oo^n,  2.0OO00OO 

P  r 

216^17, 

ft  .  ?6*Oftfti)  f 

p. ftoooono 

-2,-2  9.  ?-2£U— (1.5  2  n  o  0  4 Iv. 

2.0000000  — 

2.  ? 5  ■  2 ft  5  ■  (1.230°3ft3.  2.0000000 

2. 

?W5J  ft* 

_r;,2b8ftf)ftft. 

ft.ftOftoftnn 

2. 

2ft,  2?l  , 

0.45251*2, 

?.00onooo 

2,  25,269,  0*162*637,  2,0000000 

ft  , 

?0, *1 *, 

2.nooftooo, 

0,000 ft OftO 

2. 

?9,P?P, 

0.3ft«2316, 

2.0000000 

2,  25,?85,  0.0907592,  2.0000000 

0, 

I°.*1«, 

p  t  n  o  o  o  o  o  o  ,• 

ft .  0(»ono n o 

2. 

-J>4.  2?  3.— tt.il  *  1  “  «  ft-. 

-2*0000000  — 

_ 2,  2Sa?P6f  n. Q?n66«s._0-391 6931- 

ft  r 

lft^S! •, 

2.  Oftftftftnri , 

ft. OftOoonn 

2, 

0.20805*3, 

2.0000000 

?,  ?0,?87,  0.958S819.  2,0ftft000ft 

p. 

*^.^1  *, 

•p.ftOOOOOft  , 

ft,  oooooon 

?, 

?4 , ??5 , 

0  4  7°0717, 

?.  01*00000 

?,  24,2R»,  (1. *4(10453,  2.0(100000 

o. 

P,91  *, 

2. ftooooon , 

ft . OftOOAOO 

-2-,- 

-2Q*?26*. 

_ C_  aJB54-._2*  0000000 _ 

?.  ?4  .  ?  ft  ft  .  n.H2?4(T,7.  2.0000000- 

0.- 

_  1.-51*. 

-£.oooonoo. 

P.OOOOPftO 

2. 

2«>,  227, 

0.00370*5, 

n.*432aft5 

2,  ?4 , ?ftfl, .  0.75f 3??  ,  2.0000000 

«, 

ft  ,  9  J  A  f 

O.ponnnoo, 

o .  oooftono 

2, 

2ft.  ?2ft , 

0.97570^9. 

2. 0000000 

2,  ?0,?«M,  fl.h«*2S35,  2.0000000 

n. 

ft, 917, 

ft . ooftonoft  r 

? ,onooooo 

?T  ?*-??4.  (i-4n?*?33.  2, o o o o o o o — 

?.  24.2ft?,  11. *1*10, .4.  2.0000000 

0. 

0,51*, 

0.0000001), 

?. ooonnoo 

?, 

?*  ,  230  , 

0,83003*7, 

2.O0OOO00 

2,  20  f  293 ,  0 .5500603,  2.00000P0 

p. 

ft  ,  Xftft  , 

(  .ftOftOO  ft()  , 

P,onoftooo 

2. 

28,231  , 

o.77ioonn. 

2.0000000 

2,  ?ftf?94,'’  ft.48197’*4-,  2.C00ftft00 

p , 

o ,  «n, 

ft^onnnonr,  t 

2 :  #  0  ft  0  ft  ft  O  0 

-2.-25,  2J  2.  D.1CJ363  4, 

..2*0000000 _ 

?r  24. ?Q5 .  0.413*070.  2.0000000 

0. 

0.  B9 , 

o*opopooo. 

p. OOOOOOO 

2, 

2*. 233, 

0.f'527«ft, 

2,0000000 

2,  2«,2ft*,  0.345*0*4 ,  2.0000000 

2. 

?6,2?4, 

0.5*7 1**3, 

2.0OO0«OO 

2,  34,247,  0.2777138,  2.0000000 

-  2... 

-28,235. 

_n_ 009101 6, 

.  2.0000000 _ 

2.  24.208.  0.204*252.  2.0000000 

?, 

28,23**, 

o.ft  ft  1(1100, 

2,0000000 

2.  24,2«o,  0.14153*6,  2.0000000 

2, 

2ft. 237, 

ft. 3*202*4, 

2.00O0O0O 

2,  24,300,  0,0734521',  2.0000000 

_2,_2£,.235,_(’.2«J(  *.4  18, 

_2*  OOoooon _ 

2  r  2Q.301.  ft. 005^659.  0  07*79b9_ 

2. 

2ft. 23ft, 

0.22*7501, 

2.0000000 

2,  23,302,  0.43727*7,  2.0000000 

p,  2*,p*«o,  0.15**666,  p.n<iooooo 
.  2*-2fl,2"U—fl*0905Bn0,_2*00O00O0— 
2,  ?*f2«2,  n.n??ao^6,  0.33036*0 
?,  n.aSdflOfc?,  2.finoonnn 

-  2*_27*,Z<UU*JU.fl5632ft OOOPOPO — _ 
?,  27  ,  >4  9 ,  n.81ft?3?5,  2 . ooonooo 

2,  27,?<*6,  n.?90iasn,  2. ftOftOftOO 
_2*_2I*2a7__IL*68205B3„_  2*0000000 — . 
2,  27, 2«*,  0.b13<57l7,  2 .  onftnoon 
?,  ?7,  >uo,  0.c0S*b51,  P.opOftnnn 

_  2f_  .2?^?5Qf_  Q.rim?9*_S^_2-flflP00Q0 _ _ 

2,  2*>,25l,  0.4ft97iift,  p.onnonoo 
2 r  ?7,?52,  0.3*162*2.  ?. nnnnnnn 
. .  2*_27*.25i,— O^J-353 6 7-,. .2*  0000000 
2,  ?7,2So,  0. 205*501,  2.nonoonft 
2,  ?7  ,  ?59,  .o.  1 3T3b3<i ,  ?.f0ft0nnp 

-  2*— 2.?.*25iu _ 0*J7i-9L27fc<L*_i»  00  00000 — 

> ,  27, >57,  o.0011««2,  o.ftl7a8** 

2,  ?b, >6*,  0,.9?3JP3br  2.0000000 

.  2*  _?b#?59*_.£>*  865P150., 

2,  ?6,2bl>,  0.79*,<»;>*a, 


2,  23,303,  0.^919?1,  2.0000000 
-£,--23*.  3(LA^iUjmXfl  3  bj — 2*  0  0  0  0  0  0  C_ 

?,  23 , 3{»9 ,  0.7-«30170,  P.0000000 
2,  23,30*.,  o.bb«o303',  2.0000000, 
2.  23.307,  n  -  996*437  .  2.00OQO0Q_ 


2,  23,300, 
■2,  23,300, 
Z.  23,31 Q, 


0.5287571,  2.0000000 
0.4606705,  2.0000000 
P-3Q25B38.  2-OOftPO{K\_ 
2,  23,311,  0.3244953,  2.OO0O000 
2,  23,312,  0.2560007,  2.POOOOOP 

2^  e3,3l«;  0. 1202354,  2.noonopn 
2,  23,315-,  0.0S?1488,  0.7658O97 

2,  22.317,  0.915«7*6,  ?;0000OO0 
2,  22,31  * ,  0.0070070,  ?.r»ooclnoo 
_ 2.  ?2.3io.  f)r77o^nna. _ 2* 0  O.o r. o o o_ 


2*0000000 
2  .  AOOOAOO 

2 _ ooonono 


t  ft, >61  ,  0.72**418, 

?b*  262* — (J*b60  7_552»  — 2»  ft ft POO 00  - 
<; .  ?#rt«Oftooo 

>**,>*.«,  ft.5>*5*19,  2.0ft*ftftnn 

>«..***.  r.  »i^i.>iq<i  .  >  ftftouftnfi 


2,  22.320,  0.7117138,  2.0000000 
?,  ??,?2i,  o.fcaJfc?*2,  l.nononnn 
_ !.  '-57554*5.  2.PP*P"’’‘>_ 

?,  ?2,3?3,  ,5«7q539,  J.«0*p««n 

?,  •.«743*5«,  ?.n(i««non 

?.  ??.7?4.  r.  .ryitrn*.  ?.  nomnn'i 
?,  o. 3r>s» q?i .  ?.B»n,ino* 

2,  22.72T,  n.2>5t "«S.  2.ftn«on*r> 
2,  22.22*.  0.»*T1M«4,  ?.00fl*PP* 


Figure  1.  Scribe  Boundary  Input  Data  for  Gyrotron  Gun  (1/^mm  *€SH)  (Cent.) 
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Figure  7.  Axial  Magnetic  Field  -  Magnet  Two 
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D.  Results 


As  mentioned  earlier,  the  simulations  reported  here  were  conducted 
with  the  goal  of  determining  consistent  sets  of  gun  parameters  for  the 
successful  operation  of  the  high-power  gyrotron  now  under  design  at  NRL  for 
air  breakdown  experiments.  This  was  not  to  be  a  gun  design  exercise,  and 
no  changes  in  the  electrode  shapes  were  made.  In  addition,  extensive 
simulations  to  determine  parameter  sensitivity  around  any  one  particular 
operating  point  were  not  attempted. 


'•Vithin  the  context  of  the  parameter  constraints  listed  in  Section 
III.B.  above,  the  simulation  results  are  summarized  in  Table  2.  The 
notation  for  the  table  headings  is  as  follows: 


SCRRUN 
V 

VMOD 
BZA-K 

BZA-F 

ALPHA 
R  G-C 
BETA  SP. 

The  table  entries  are  listed  in  the  following  order:  l-decreasing 

accelerating  voltage  (anode  voltage);  for  each  accelerating  voltage, 
2-decreasing  magnetic  field  at  cathode;  3-decreasing  mod  or  intermediate 
anode  voltage.  Based  on  the  criteria  listed  above,  several  parameter  sets 
have  been  selected  as  possible  gun  ooeratinq  point's.  ,  These  are  indicated 
in  Table  2  by  asterisks  (*).  Those  simulations  using  the  1/4mm  mesh  are 
indicated  by  UOTE  5. 

The  ciun  perveance  is  fixed  in  these  simulations  to  model 
tenperature-limited  emission  at  10A.  However,  two  runs  (11,12)  were  made 
at  reduced  oerv^ance  to  simulate  emission  at  5.9S  and;7.9^A,  respect ivel\ . 


SCRIBE  run  identification  no. 

anode  voltage  (kV) 

intermediate  anode  voltage  (kV) 
on-axis  axial  magnetic  field  at 
cathode  axial  locaton  (kG) 

on-axis  axial  magnetic  field  in 
uniform  field  region  (kG) 

transverse/axial  velocity  ratio 

average  guiding  center  radius  (cm) 

spread  in  relativistic  factor  beta  (°o) 


© 


* 


SCRUM. 

V 

VI10D 

BZA-k 

9 

71 

35 

3.73 

10 

71 

3S 

2.9 

96 

70 

34 

2. SI 

95 

70 

35 

.  2.76 

93 

70 

37.1 

2.76 

*  SI 

70 

37 

2.5 

*  SO 

70 

35 

2.5 

*  6 

70 

33.05 

2.37 

*  2 

70 

33 

2.37 

*  7S  ' 

70 

35 

2.3 

*  79 

70 

32.9 

2.3 

*  S3 

67 

35 

2.5 

*  S4 

67 

35 

2.45 

*  55 

67 

35 

2.4 

11 

65 

35 

2.9 

52 

65 

34 

2.5 

14 

60 

35 

3.17 

13 

60 

35 

2.9 

12 

60 

35 

2.9 

*  60 

60 

35 

2.37 

*  63 

60 

Vx 

2.37 

*  61  • 

60 

33 

2.37 

*  '62 

60 

32  • 

2.37 

*  24 

60 

'  30.5 

2.37 

20 

60 

30 

2 .37 

44 

60 

30 

2.37 

*  77 

60  . 

30.5 

2.3 

*  75 

60 

30 

2.3 

*  76 

60 

30 

2.3 

B2A-F  ALPHA  R  G-C  BETA  SP. 


33.66 

.473 

.241 

1.6 

33.59 

.762 

.205 

7.6 

33.34 

.539 

.204 

4.4 

33.57 

.S75 

.197 

4.5 

33.57 

.532 

.199 

5 

33.55 

1 .39 

.1S3, 

3.31 

33.55 

1  .13 

.191 

3.21 

33.54 

1.24 

.171 

3.9 

33.54 

1.24 

.17 

4.2 

33.53 

2.31 

.167 

2.55 

33.53 

1.51 

.179 

5.41. 

33.55 

1.21 

.196 

4.09 

33.54 

1.39 

.193 

4 .45 

33.54 

1 .65 

.189 

5.29 

33.59 

.77 

.174' 

9.6 

33.55 

1.1 

.191 

3.25 

33.61 

.791 

.214 

5.1 

33.59 

.S49 

.2 

12.4 

33.59 

.655 

.2 

11 .5 

33.54 

2.16 

.195 

4.45 

33.54 

1.31 

.186 

4.23 

33.54 

1  .64 

.196 

6.3 

33.54 

1 .39 

.199 

6 

33.54 

1.13 

.174 

4.2 

33.54 

1 .06 

.IS 

4.3 

33,54 

1  .OS 

.196 

5.9 

33.53 

1.39 

.179 

4.72 

33.53 

1  .27 

.14 

5.1 

33.53 

1.3 

.147 

6. S3 

NOTE 


1 


2 

2 


2 

2 

5 

5 

5 

3 

2 


4 


5 

2 

2 

3 


*  -  POSSIBLE  OPERATING  POINT 
Notes: 

1  -  Slight  Origin  Shift  for  Maanetlc  Field 

2  -  Step  Size  Halved 

3  -  4.95A  ' 

4  -  7.94A 

5  -  Iwice  Original  Resolution  (1/4mrt) 


Tahle  2.  Simulation  Results 
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v/l  U1  U’ 


9 


9 


4 


% 


4 


The  qun  performance  does  not  appear  to  be  extremely  sensitive  to  the  rela¬ 
tive  positions  of  the  qun  and  coil  systems,  as  evidenced  by  a  comparison  of 
runs  95  and  96.  The  coil  system  was  shifted . axially  bv  approximately  0.9cm 
between  runs  95  and  96.  Hith  the  exception  of  run  96,  the  magnetic  field 
at  the  z  coordinate  of  the  emission  surface  scales  as  Bz(kG)  -  -0.9558  1(A) 
+  4.5091,  where  I  is  the  trim  coil  current  and  with  the  main  coil  current 
fixed  at  45.41A. 

,  Representative  trajectory  and  trajectory  with  equipotential  plots 
from  these  simulations,  including  both  the  1/2  and  1/4mm  meshes,  are  shown 
in  Figures  S  -  15.  These  plots  are  presented  in  the  same  order  as  listed 
in  Table  2.  Since  azimuthal  symmetry,  is  assumed  for  these  simulations, 
only  the  first  quadrant  of  the  electron  qun  is  shown.  The  positions  of  the 
cathode,  intermediate  anode,  and  anode  are  indicated  in  Figure  8.  Note 
that  the  "R"'  axis  scale  is  approximately  5.6  times  that  for  the  "Z"  axis 
scale  for  enhanced  visual  cxarity.  The  scales  for  both  axes  are  in  mesh 
units.  Six  trajectories  are  followed  with  the  1/2mm  mesh,  and  ten  are  fol¬ 
lowed  with  the  1 /4mm  mesh.  The  electron  emission  takes  place  approximately 
1.5  cells  from  the  actual  emission  surface.  The  approximate  center  of  the 
actual  emission  surface  is  located  at  (0.6cm,  1.1cm). 

Although  simulations  of  the  full  parameter  matrix  were  not  made, 
enouqh  data  exists  to  demonstrate  several  trends.  The  most  important 
parameters  for  successful  gvrotron  operation  are  alpha  and  the  average 
guiding  center  radius.  '  The  variation  of  alpha  with  the  intermediate  anode 
voltaoe  for  constant  anode  voltage  and  magnetic  field  is  shown  in  Figure 
16.  In  this  case,  V=60kV ,  the  magnetic  field  at  the  cathode  is  2.37  kG, 
and  the  final  magnetic  field  is  33.54  kG.  Although  typical  power  supply 
construction  does  not  allow  a  continuous  variation  of  the  intermediate 
anode  voltage,  one  can  see  that  this  would  be  a  convenient  'knob'  to  turn 
to  adjust  alpha.  Under  the  same  conditions  as  in  effect  for  Figure  16, 
relatively  small  variations  in  intermediate  anode  voltage  do  pot  have  a 
significant  effect  on  the  average  guiding  center  radius.  This  is  shown  in 
Figure  17. 
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Figure  6.  Trajectories  and  Equipotentials:  Run  96 
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Trajectories:  Run  90 


27 


filROTRON 


5GRRUN.S3R 


Figure  10.  Trajectories:  Run  S3 


<wm  »/nsixu  v 
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Figure  14.  Trajectories  and  Equipotentiais: .  Kiin  4S 
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A  convenient  experimental  parameter,  however,  is  the  magnetic 
f;'eld  at  the  c, .  :node.  Figure  18  shows  the  scaling  of  alpha  with  cathode 
magnetic  fiald  for  two  different  anode/intermediate  anode  voltage  settings, 
67kV/35kV  aid  60kV/35kV.  The  spread  in  beta  (parallel  velocity/speed  of 
light;  f  .  tne  former  case  is  shown  in  Figure  19. 

In  summary,  a  number  of  simulations  of  gvrotron  electron  gun 
performance  under  a  given  set  of  conditions  have  been  made.  The  results  of 
these  simulations  have  shown  that  several  possible  operating  points  do 
exist,  and  that  the  operating  conditions  vary  approximately  linearly  for 
relatively  small  changes  in  the  typical  operating  parameters. 
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Figure  19,  Velocity  Spread  Scallop  with  Catnode  Magnetic  Field 


REFERENCES 


1.  A.  0.  MacDonald,  Microwave  Breakdown  in  Gases,  Wiley,  New  York, 

1966.  ■  ‘  “  ~ 

2.  This  technique  is  also  used  to  investigate  plasma-enhanced  thermal 
,  coupling  to  metal  surfaces  by  T.  3.  Wieting,  3.  L.  DeRosa  and  M.  E. 

Read  in  a  parallel  program. 

3.  C.  L.  Yee  and  A.  W.  Mi,  "Microwave  Energy  Deposition,  Breakdown  and 
Heating  of  Nitrogen  and  Air,"  NRL  Memo  Report,  No.  4617,  1981. 

4.  W.  W.  Salisbury  and  W.  Flynt,  "Microwave  Energy  Coupling  Through 
Re-entry  Plasma,"  Varo,  Inc.,,  RADC-TR-64-575,  December  1964. 

5.  W.  Beust  and  W.  L.  Ford,  Microwave  3ournal,  M10),  p.  91,  1961. 

6.  Y.  P.  Raizer,  Sov  Phy  3ETP,  34,  p.  114,  1972. 

7.  W.  E.  Scharfman,  W.  C.  Taylor,  and  T.  Morita,  IEEE-Trans.  Antennas 
Propag.,  12,  p.  709,  1964. 


S.  M.  E.  Read,  R.  M.  Gilgenback,  R.  Lucey,  K.  R.  Chu,  A.  T.  Urobot, 

,  "Spatial  and  Temporal  Coherence  of  a  35  GHZ  Gyromonotron  Using  the 
TEci  Circular  Mode,"  IEEE  Trans.  Microwave  Theory  Tech.,  28,  p.  875, 
1950.  . 

9.  W.  R.  Ott,  P.  Fieffe-Prevost,  and  W.'  L.  Wiese,  Appi.  Opt.  12*  p.  1615, 
1973. 

10.  C.  R.  Barber,  3.  Sci.  Instrum.,  23 ,  p.  235,  1946. 

11.  H.  R.  Griem,  Plasma  Spectrostrophy ,  McGraw-Hill,  New  York,  1964. 

12.  M.  Imami  and  'V.  Brost,  3.  Chem  Phys.,  61 ,  p.  1115,  1974. 

13.  R.  M.  St.  3ohn,  F.  L.  Miller,  and  C.  C.  Lin,  Phys.  Rev.  A,  134, 

p.  A888S,  1964.  - 

14.  R.  Nichols,  Ann  Geophvs.,  2£,  p.  144,  1964. 

15.  w.  L.  Wiese,  M.  Smith,  ang  B.  M.  Clenn,  "Atomic  Transition 

Probabilites,"  Vol.  1,  National  Bureau  of  Standards  Reference  Data 

Systems  (NSRDS-N9S-4) ,  1966. 

16.  P.  Millet,  Y.  Salamero,  H.  Brunet,  3.  H.  Galy  and  D.  Blanc,  3.  Chem. 

Pnys.,  58,  p..  5-5  39  ,  1973. 


3? 


REFERENCES  (Cont.) 


VV.  M.  Bollen,  R.  H.  3ackson,  B.  Goplen,  R.  3.  Barker,  R.  E.  Clark, 
3.  McDonald,  C.  Sedlak,  "Computational  Capabilities  For  Microwave 
Applications,"  MRC  Report  MRC/WDC-R-OSO ,  March.  1934. 

M.  E.  Read,  private  communication. 

S.  3.  Sackett,  "EFFI  -  A  Code  for  Calculating  the  Electromagnetic 
Field,  Force,  and  Inductance  in  Coil  Systems  of  Arbitrary  Geometry," 
Lawrence  Livermore  Laboratory  Report,  UCRL-5240^,  March  1978. 


FILMED 

10-85  ' 

DTIC 


